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Abstract 

We study a model of 2D QFT with boundary interaction, in which 
two-component massless Bose field is constrained to a circle at the 
boundary. We argue that this model is integrable at two values 
of the topological angle, 6 = and 9 = ir. For 8 = we propose 
exact partition function in terms of solutions of ordinary linear 
differential equation. The circular brane model is equivalent to 
the model of quantum Brownian dynamics commonly used in de- 
scribing the Coulomb charging in quantum dots, in the limit of 
small dimensionless resistance go of the tunneling contact. Our 
proposal translates to partition function of this model at integer 
charge. 
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1 The model 



In this note we address a model of 2D Euclidean field theory with boundary 
interaction which we call the circular brane model. The model contains two- 
component Bose field ~X.(z,z) = (X 1 (z, z), X 2 (z, z)) defined on a disk of a 
radius R with usual complex coordinates (z,z), \z\ < R. Its bulk dynamics 
is that of a free massless field, as described by the action 



Aiisk = - f d 2 zd z x-d- z x (i) 

l\z\<R 



IT 



(d 2 z = dxdy), but it obeys nonlinear boundary condition 

X 2 B = x , (2) 

where Xg stands for the values of X at the boundary \z\ = R, and go is 
a coupling constant. Due to the nonlinear boundary condition, the theory 
needs renormalization. It has to be equipped with the ultraviolet (UV) 
cutoff, and consistent removal of the UV divergences requires that the bare 
coupling constant go be given a dependence of the cutoff momentum A, 
according to the Renormalization Group (RG) flow equation 

A^ = -2 5o 2 -4<7o 3 + --- ■ (3) 

The leading two terms of the beta function written down in (J3J) were com- 
puted in m and |2j, and indeed agree with the more general calculations 
in [3] and [I] . Eq. (j3J) shows that the theory is asymptotically free - at short 
distances the effective size of the circle (|2| becomes large. It is believed that 
the theory develops a physical energy scale E* , and at the distances 3> 1/E* 
from the boundary the effect of the boundary is that of the fixed boundary 
condition = (the boundary @ "flows" to the Dirichlet boundary). 
According to the equation Q 

E* ~ A g^e"^ . (4) 

As usual, the effect of the boundary can be described in terms of the 
boundary state 0- The boundary state \B) is a special vector in the 
space of states of radial quantization (in our case the space of states of two- 
component free massless Bose field on a circle of the circumference 2nR). 
This state incorporates all information about the boundary conditions; any 
correlation function can be written in terms of the Euclidean-time Heisen- 
berg operators sandwiched between | B ) and the radial-quantization vacuum 



1 



I ). In particular, the overlap ( | B ) coincides with the disk partition func- 
tion, up to the factor i? c / 6 . In our case | 0} is the Fock vacuum of the free 
massless field X(z, z) with the zero-mode momentum P = 0. More gen- 
erally, we will be interested in the overlap amplitude ( P | B } , where | P ) 
is the Fock vacuum with the zero-mode momentum P; of course, this am- 
plitude relates to the (unnormalized) one point correlation function of the 
exponential field inserted at the center of the disk, 

(P|£) =j R-/6+P 2 /2 (e iP-X (0)0))disk _ (5) 

Here c = 2 is the central charge of the bulk theory. It will be convenient 
to transform the spherical brane model to the geometry of semi-infinite 
cylinder, using the standard exponential map: 

Z /R = e u / R , z/R = e a/R , (6) 

where u = a + ir, u = a — ir, with a ranging from — oo to 0, and r being 
a periodic coordinate with the period 2nR. In what follows we will write 
(with some abuse of notations) X(r, a) for the field X in expressed in 
terms of these cylindrical coordinates. Then, introducing the shifted field 
variable 

Y(t,<t)=X(t,<t) + ^P , (7) 



one can relate the the amplitude (P\B) to the partition function 



i. 



Z(R,H) = Km |2ess J VY e ~ Ac y l[Y] j , (8) 

where the integration field variable Y(r, a) is still subject to the circular 
brane constraint 



and 



Ayi[Y] = i- J daj&r (d a Y ■ d a Y + d T Y ■ d T Y) (10) 



2 <b dr H • Y b ■ 



lr The factor 2 in @ is added partly to eliminate the spurious boundary degeneracy 
of the auxiliary Neumann boundary at a = —L in il l lit , and partly to take into account 
the infrared boundary entropy of the real boundary at a = 0; with this normalization the 
R — *• oo behavior of Z(R, H) is ~ e~ 27rR£ , with the coefficient 1. 



2 



Here and below § dr stands for ^ R dr. The "external field" H relates to 
P in © as 

and the argument H in the r.h.s. of (JHJ) is \/H 2 . Thus 

(P\B) = ±Z(R,^) , (12) 

where the factor is the boundary degeneracy of the infrared (Dirichlet) 
fixed point. We note that by the nature of the boundary condition (jUJ) the 
partition function Z(R, H) is expected to be an entire function of H, and its 
analytic continuation to pure imaginary H is unambiguous. In the form (|10[) 
the model has interpretation as 1+1 QFT on a half-line, at the temperature 
T = ^r- At R — > oo the partition function (JSJ) develops standard linear 
asymptotic 

log Z(R,H) -> -2ttR £(H) as R -» oo , (13) 

where £ (IT) is the boundary energy. This quantity depends on H and <7o> 
the latter dependence being through the energy scale E* , Eq.(j3J). We found 
it convenient to fix overall normalization of E* relating it to the zero tem- 
perature susceptibility, 



1 1 d 2 £ 



E* 2 dH 2 



(14) 



The circular brane model shows many similarities with the 0(3) sigma 
model (or n-field), the asymptotic freedom © being just one of them [7]. 
Another similarity is apparent when one observes that the field configura- 
tions Y(t, a) in (JHJ) can be separated into the topological classes character- 
ized by the integer- valued winding number 

w= £ f dr Yb a9tYb ' ( is ) 

just like the sigma-model field configurations n fall into their topological 
classes characterized by the mapping degree S 2 — ► S 2 . Moreover, in the 
sigma-model the instanton configurations minimizing the action in each 
topological sector are found explicitly in terms of rational functions [S]; 
very similar expressions exist for the instanton configurations in the circu- 
lar brane model |Hlll()j. The calculations of the instanton determinants jllj 
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also exhibit much similarity |121ll3llTl] . As usual, existence of the integer- 
valued topological charge allows one to introduce another parameter, the 
topological angle 9, 

Z e = Y, z(w) ei9w , ( 16 ) 

where the partition functions Z^ w ' are the functional integrals (jSJ) taken over 
the field configurations in the topological class w only. The full partition 
function Z, as defined by Eq.(JSJ), involves all topological sectors with equal 
weights, i.e. Z = Zq = q. It is this 9 = partition function which we had in 
mind in defining the boundary energy £{H) and the energy scale E* through 

One of the reasons we mention this analogy here is because the 0(3) 
sigma model is known to be integrable at two special values of the topological 
angle, 9 = and 9 = ir |15 |I16[IT?| I18|. Therefore it is natural to expect that 
the circular brane model is also integrable at these two special values of 9 2 . 
We argue in Appendix that for these values of 9 the model admits a number 
of nontrivial higher-spin local integrals of motion. In fact, our argument 
there goes for the whole class of the u O(N) spherical brane models", with 
TV-component field X, subject to the spherical brane boundary condition 
(j2|. Moreover, in these models it is possible to describe the whole set of the 
local integrals of motion l2fc-i, k = 1, 2, 3 . . . in some details. 

The aim of this note is to present a proposal for the exact form of the 
partition function Q at the integrable point 9 = in terms of solutions 
of certain ordinary differential equation. We will present this in the next 
section, but let us first make some remark and introduce suitable notations. 
As was mentioned above, renormalization trades the bare coupling constant 
go for the RG invariant scale (|4"|). and the partition function actually de- 
pends on the dimensionless combination RE* . In fact, in the circular brane 
model (as well as in the 0(3) sigma model) the last statement is not entirely 
true. As was observed in |12j . the functional integral (JSJ) has a specific non- 
perturbative divergence due to the small-size instantons, which cannot be 
absorbed into the renormalization of the coupling constant go- Due to this 
effect, the partition function (JSJ) is expected to have the form 

Z e = e LKCOsd Z* 9 ( K ,h) , (17) 

2 It is worth mentioning here that the case 8 = tt is of special interest. In this case 
the circular brane model flows to a nontrivial conformal boundary condition which is 
believed to be related to the infrared fixed point of two-channel spin-| Kondo model in 
its interacting sector |19| . 
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where 



L = A log (B (g ) A/ E*) , (18) 

and 

K = 2irRE* , h = 2ttRH . (19) 

Here A is a constant, which is universal (actually A = 2), but B(go) is not - 
it depends on the details of UV cutoff. The factor Zq depends on the RG in- 
variant parameters k, h and 9 only. Similar nonperturbative small-instanton 
divergence is well known to be present in the 0(3) sigma model where 
it is blamed, for instance, for the violation of the normal scaling of the topo- 
logical susceptibility j2Hj. It is important to note that the small instanton 
divergence does not affect the /i-dependence of the free energy (interaction of 
small instantons with the external field H is negligible) , therefore the energy 
scale E* defined through ()13j) and (|14|) enjoys the normal RG behavior 
In what follows we discuss the universal factor Z*(k, h) in the full partition 
function at 6 = 0, 

Z(R,H) = e KL Z*(K,h), (20) 

which carries all interesting dependence on H and R. Since any factor of 
the form e const K can be absorbed into the redefinition of B(g$) in (|18j). we 
fix the ambiguity by assuming that Z*(k, 0) — > 1 as k — > oo. With this 
convention, the boundary energy can be written as 

S(H) = LE* + £*(H) , (21) 

where £*(0) = 0. 

2 The partition function 

It was discovered some years ago in Ref. |21j that in the case of the boundary 
flow in the minimal CFT with integrable boundary perturbation |22| . the 
overlap analogous to © can be related exactly to the eigenvalue problem 
of certain ordinary differential operator. Despite the fact that this relation 
was proven [2^] , and similar relations were found in other integrable models 
of CFT with non-conformal boundary interactions |24|I25| . its deeper reason 
remains mystery to us. Nonetheless, it looks natural to assume that it might 
be a general phenomenon, and try to identify associated differential operator 
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for the circular brane model. Below we simply present our proposals for 
such differential operator. The motivations came from detailed studies of 
somewhat more general, but still integrable "brane" model (we call it the 
"paperclip brane model"); it generalizes the circular model in a way similar 
to the "sausage" deformation of the 0(3) sigma model [SHj. We will report 
these studies elsewhere [27] . 

Consider the following ordinary differential equation 

{ -d 2 + K 2 exp (e v ) + /i V } 9(y) = . (22) 

Let ip-( v ) be the solution which decays at large negative v, 

ip-(v)^e KV as v^-oo. (23) 
Also, let 4>+{v) be the solution decaying at large positive v, 

ip+(v) — > exp < — — + k Ei^— J > as v — > +oo , (24) 

where Ei(x) is the integral exponent function, Ei(x) = P.V. f^° x 4- e_i - The 
asymptotics (j2*3*|) and (|2*I|) specify the solutions tp-(v) and ip+(v) uniquely, 
including their normalizations. Then 

l~WTW w[ * + '*- ] ' (251 

where W is the Wronskian ip + d v ip- — ip- d v ip + and = 0.577216 ... is 
Euler's constant. Thus determination of Z* reduces to the problem of the 
ordinary differential equation. 



3 Low temperature expansion 

At large R (low temperatures) the equation (f2*2*|) can be studied using WKB 
expansion. The leading WKB approximation yields explicit expression for 
the boundary energy (|21jl. 

/•OO Af , \ 

£*{H) = -E* J j (y/tf + HH/E*) 2 -^ 2 } . (26) 
Note that the integral admits the expansion in powers of H 2 , 

oo 

£*(H) = E* (#) 2 \ (27) 

k=i 
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with 



This expansion should be understood in terms of the expansion (|66|) . with 
the powers H 2k representing the highest order terms h 2k of the polynomials 
I^k—ii the vacuum eigenvalues of the local IM Igfc—i- Note that the above 
expression for the coefficients Cj. is consistent with (|71|l. The next order 
of the WKB expansion provides leading correction to the linear asymptotic 

CEP, 

logZ* = -2kRS* -JLs* + 0(R- 3 ) , (28) 
where, by direct calculation, 



1 d_f 

8H dH\ ~dH 

Again, expanding l|28[) in powers of H we have 

oo 

logZ* ~ -J2 C k (2nE*) 1 - 2k I 2k _ x , (30) 



k=l 



with 



in perfect agreement with Eas. (|66|) . (|69j) . Further terms in (|31|) can be 
obtained by computing yet higher orders of the WKB expansion. This way 
one finds 



R V 12 



i? 5 V 4 100 5040 
again in agreement with (|68|), 
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4 High temperature expansion 



It is also possible to develop short distance (high temperature) expansion of 
the partition function (|25|). For small n the solutions ip-{v) and ^+{v) can 
be evaluated using appropriate versions of perturbation theory. To make 
the resulting formulae compact we shall describe here these expansions in 
the case h = only. The solution (v) can be written as 

oo 

t/>-(v)=e KV {l + J2 a n(K)e nv }, (33) 

n=l 

where the coefficients a n (K) = 0(k 2 ) and admit systematic expansions in 
powers of k. The k — > expansion of the solution ip+(v) is more subtle. In- 
stead of powers of k, it expands in powers of the "running coupling constant" 
g = g(n) defined by the equation 

k = g~ l e~2i . (34) 

After change of variables 

v = g x - log(#) (35) 

the equation (I22H becomes 

{ -dl + e x + 5U(x)}^ = 0, (36) 

where 

ffZ7(x) = exp ( j - e x . (37) 

For |x| ~ 1 the term 5U(x) is 0(g); in this domain solutions of (|3fij) ad- 
mit systematic expansion in powers of g. Expansion of ip+ is obtained by 
iterations of (|36j) . starting with ip^/ = C(g) K (2e x ^ 2 ) , where K is the 
Macdonald function and the normalization constant C(g) must be adjusted 
to match the asymptotic form (|24|). It turns out that both expansions, of 
ip— in the powers of k, and of V'+ in the powers of g, have common domain of 
validity at v ~ — log(g), and can be used there for evaluation of the Wron- 
skian in (|25jl . As the result, the following form of the small k expansion of 
the partition function (|2T)|) emerges 

K 00 

Z(R, 0) ~ e 2K l°g(2-B(90)i?A) ^_ y K n ( 3g ) 
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where the coefficients z n (g) are power series in g. The term n = in (|38jl 
should be interpreted as the perturbative contribution to the partition func- 
tion, the term n = 1 corresponds to the one-instanton contribution, and 
the higher nonperturbative terms n = 2, 3 . . . presumably describe the con- 
tributions of multiple instanton and anti-instanton configurations. Explicit 
computation along the lines described above yields 

zo(g) = 1- (1 + je) 9 + 0(g 2 ), (39) 
Zl {g) = \og{2)- 2 + 2, lE - 2 lE g - 2{2 + lE ) g 2 + 0{g i ) . 



5 Dissipative quantum rotator 

The circular brane model has useful interpretation in terms of Brownian 
dynamics of quantum rotator. It was noticed a while ago in Ref. |2E] that 
the free massless bulk dynamics (fTU)) is equivalent to the Caldeira-Leggett 
model of quantum thermostat |29j . Upon fixing the boundary values Yb(t) 
and integrating out the bulk part of the field Y(r, a), Eq.(JHJ) reduces to 



Z(R,H) = J Vrj e - Adiss[v] , (40) 

with 

2H f 

Aiissh] = ■= <p dr cos(ry) + (41) 

\/9o J 



1 f / sin 
dr <Z> dr 



2 ( v(t)-v(t') 
2 



8n 2 &go J J sin 2 (^) 

where t](t) is the angular field defined through H-Yb(t) = H/ '■ s /g~o cos r/(r). 
The model similar 3 to Q41JI was introduced in (see e.g. for a review) 
as an effective field theory describing tunneling of quasiparticles between 
superconductors. More recently it was subject of much interest in studying 
of the phenomenon of Coulomb blockade in quantum dots jS2j , in the regime 
where the dimensionless resistance of the tunneling contact go is small |33( 
I12l l2i rn?lll4| . In this context, the effective action in fact contains another 
term, 

Acb = A diss + j^- j> dr rfr , (42) 



3 The model of |3U| has the potential term cos(2r7) instead of cos(?7) in 14111 . 
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where Ec has the meaning of the charging energy of the dot in the absence 
of the tunneling, E c = e 2 /2C. Of course, r is the Matsubara time, and 
2ttR = 1/T, the inverse temperature. When T <C Eq this term just provides 
explicit UV cutoff, with the cutoff energy 

A = Ec/go ■ (43) 

With the cutoff procedure thus specified, the "small-instanton" factor e nL 
in (|20[) becomes unambiguous. By comparing (|38[) . (|39|) with direct one- 
instanton calculation in 1)42(1 |12U13llT4*| one finds for small go 

L = ^ + 5 1og( 50 ) + O(l) • (44) 

Eas.(|20 j) .(|25 |) then provide the partition function of the model (|4*2*)) at 6 = 0, 
in the limit of small go . One only has to relate the energy scale E* in (|19|) to 
the parameters go and Ec in the above action (fl2*|) . When Ec/go ^> H E* 
and go is small, the boundary energy £{H) (|21|) can be computed directly 
from (|42[). using the standard perturbation theory in go- By developing 
two-loop order in this perturbative expansion and comparing it with the 
H/E* ^> 1 behavior of the integral in (|26|) . we found 

Ec r , 3tt 



E* = 7T~5~2 e 230 



l-—go + 0{gl)}. (45) 



Finally, let us make a remark on the topological susceptibility in the 
model ()42jl . As was already mentioned, the topological susceptibility 

E h = -iw^«)\ m («> 

in the circular brane model suffers from the small-instanton divergence, 
which leads to the non- universal factor e KL cose in (|17|) . On the other hand, 
in applications to the Coulomb blockade, the topological susceptibility in 
(|42j) has direct physical meaning - it coincides with measurable capacitive 
energy of the quantum dot. Explicit cutoff in the action (|42j) makes the 
small-instanton factor unambiguous. According to (|17|) . at zero tempera- 
ture and zero H, E c has the form 2-7T 2 (L + const) E*, where the constant 
comes from the ^-dependence of the second factor in (|17|). Since our pro- 
posal is strictly limited to the case 9 = 0, it does not allow for determination 
of this constant. However, at small go the dominating contribution comes 
from the term LE* , and therefore in view of (|44j) and (|45l) . the charging 
energy at zero temperature is estimated as 

Ec 

go 1 90 



E* c \ T =o = — e 2 «o 



l [l + 5i og5Q + o(i)] . (47) 
On > 
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Of course, the leading behavior (1/go in the brackets) here is consistent 
with the one-instanton analysis in the Gaussian approximation |12U13llTlj . 
Nontrivial prediction which follows from our proposal is the logarithmic term 
in (jl7|t . 

6 On other spherical brane models 

The circular brane model is a particular case N = 2 of the "O(JV) spherical 
brane model", the latter containing iV-component field X, subject to the 
same boundary condition (J2J), so that X# G § N ~ l 1 j. As we argue in 
Appendix, the spherical brane model is integrable for any JV > 1. Although 
the bulk part of this paper is devoted to the circular case JV = 2, here we 
would like to say few words about other cases. First, for generic JV ^ 2 
there are no topological sectors, and the small-instanton divergence is not 
present. For JV ^ 2 the overlap (JSJ) (c = JV) can be written as 

(P\B) = 2~~ e- KCo Z*( K ,h), (48) 

where k and h are still the dimensionless parameters (j!9j) . and E* is defined 
in Ea. (|14[) . The factor Z*(k, h) here is normalized as in (|2(jp. i.e. Z*(k, 0) — ► 
1 as k — ► oo, and the JV-dependent constant Co is the same as in Eqs. 

([HoT>.(|7ni). 

For generic we were not able to find differential equations which could 
possibly generate the overlap (|4Hjl . as ()22D does for N = 2. However, for two 
other cases, iV = 1 and N = 3, representations of the function Z* in (|48[) in 
terms of solutions of certain differential equations exist. Here we would like 
to describe what these differential equations are. 

For iV = 1, it is useful to note that the model is equivalent to the inter- 
acting sector of the one-channel spin-^ Kondo model. At JV = 1 the sphere 
degenerates to two points, which represent two states of the impurity spin. 
Exact equivalence is established by bosonizing the Kondo fermions (see e.g. 
Ref. [34 ). As is well known, the Kondo model is solvable by means of Bethe 
ansatz technique [3"E|l3fi| . Alternative (and equivalent) representation of the 
partition function can be given in terms of the solutions of the differential 
equation [221 

{ - d 2 v + 2vtk 2 ve v -h 2 } *(u) = . (49) 

The function Z* coincides with associated Stokes multiplier (see |25| for the 
details). 
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For N = 3, consider the differential equation 



{ - d 2 v + ttk 2 e v2/2 + /i 2 } = . (50) 
Let tp+ and ip- be its solutions uniquely determined by the asymptotics 

ip + — > exp | — ^ — ttk Erfi(u/2)j as v — > +oo , (51) 

^_ — > exp | — ^ — 7tk Erfi(— v/2) j 



,2 

as v — ► — oo . 



where Erfi(x) is the imaginary error function, Erfi(x) = 4= J^dte* . The 
suggested relation is 

Z* (k, h) = (2«VSr) -1 W [ , ^_ ] . (52) 



Like in the case N = 2, the k — > oo expansion of 1)52(1 can be obtained 
by WKB analysis of the solutions of (|51|). Thus, the leading WKB ap- 
proximation yields the boundary energy (|70j) . with the coefficients given 
exactly by (|71|) with = 3. Moreover, higher orders of the WKB expansion 
reproduce the eigenvalues l|68 jl .l|69j l perfectly. 
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7 Appendix 

7.1 The 0(N) spherical brane model 

In this appendix we present arguments in favor of integrability of the circular 
brane model at 9 = and 9 = ir. In fact, the arguments apply without much 
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modification to more general model, which contains iV-component Bose field 
X = (X 1 , • • ■ ,X N ^j, with N > 1, but otherwise is described by the same 
action (^Q) and the "spherical brane" boundary condition (j2J). Therefore 
throughout this appendix we mostly address this iV-component spherical 
brane model. Presence of arbitrary number of components allows for a 
number of checks within the 1/N expansion. 

7.2 Integrals of motion 

In the bulk, the spherical brane model is a free Bose field, and as such it 
certainly has infinite number of integrals of motion. The derivative 9X, 
where d = d u stands for the derivative over u = a + ir, is a holomorphic 
field in the bulk. Hence any local polynomial P(u) = P(<9X, <9 2 X, ■ • • ) of 
this and higher derivatives with respect to u gives rise to a local integral of 
motion 

I[P]= <f^P(d*,d 2 X,...) , (53) 

where we have assumed the coordinate a along the cylinder to be the Eu- 
clidean time. Here and below the bulk composite fields entering these poly- 
nomials are always understood in terms of the normal ordering with respect 
to the standard Wick pairing ( X a (ii, u) X 6 (u', u') ) = —5 ab log \u — u'\ cor- 
responding to a bulk free field in an infinite space. Certainly, there is a 
"left-moving" counterpart to each of the integrals (|o^|) , 

I[ j P]=/^^X,a 2 X,...) , (54) 

where d = du- The integrals (|5l^l) . ((5^1) are operators acting in the space of 
states of free iV-component massless Bose field quantized on a spatial circle 
of the circumference 2ttR, i.e. 

H = j Fp <g> Fp , (55) 

where Tp is the Fock space of the right-moving bosons with the zero-mode 
momentum P. Large set of integrals (|53jl contains many infinite subsets 
I[-P s ] (corresponding to special sequences of the polynomials P s ) of mutually 
commutative operators, 

[I[P S ],I[P S ,])=0 . (56) 
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The polynomial fields P s can always be chosen to have definite spin, and 
it is convenient to label them accordingly; below we always assume that s 
indicates the spin of P s (u). We will also use conventional notation I s _i for 
the integral I[P S ]. 

As is discussed in j37] , integrability of the field theory (JJ) in the presence 
of the boundary requires that the boundary state satisfies the equations 



(I s -I s ) \ B) = , 



(57) 



for all I s in one such subset. These equations must follow from special 
properties of the boundary conditions, which should be such that the cor- 
responding fields P s (u) and P s (u), being brought to the boundary a = 0, 
satisfy the conditions 



(P s (u)-P s (u) 



(58) 



where O s (t) are some (renormalized) local boundary fields. In the case 
of the spherical brane boundary condition ([%]). existence of the first few 
polynomials satisfying (|58|) can be demonstrated by an argument similar 
to that previously used to support integrability of the 0(N) sigma model 
[TB] . The O(N) symmetry of the boundary condition ((21) suggests that the 
polynomials P s must have this symmetry as well. We thus can look for the 
corresponding integrals I s _i in the form: 



Il = 
I 3 = 
I 5 = 

and for generic k 
I2/C-1 



dr 



2tt 
dr 
2^ 
dr 
2^ L 



ax • ax . 



(ax -ax) 2 + 63 (a 2 x-a 2 x) 



(ax • ax) 3 + 65 (a 2 x • a 2 x) (ax • ax) 



(59) 



+c 5 (a 2 x • ax) 2 + d 5 (a 3 x • a 3 x) 



j> ^ [(ax • ax) fc + 62*-! (<9 2 x • a 2 x) (ax • ax 



fc-2 



+ C2fc _!(a 2 x • ax) z (ax • ax) fe " 3 + . . . 



(60) 



Note that in writing these expressions we have fixed the normalization of the 
currents P2fc : the term (ax • d~K) k , having the highest power of X, comes 
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with the coefficient 1. In what follows we always assume this normaliza- 
tion of the integrals lafc— i- The terms omitted in (j6U|) have powers of X 
equal 2k — 4 or lower. In general case, the difference (P s (u) — P s (u) ) | Q 
is a combination of local boundary fields having appropriate symmetries. 
These include the O(N) symmetry and the anti-symmetry with respect to 
the reflection r — > — r. Note that for N = 2 the last reflection symmetry 
is violated unless 6 = and 9 = ir. Also, only the fields of right scale 
dimensions (i.e. equal or below the scale dimension of P s ) are admitted. 
Elementary counting shows that for s = 2 the only admissible field is a total 
derivative over r, of the form of the r.h.s. of Q58JI . For s = 4 there is only 
one admissible field which is not a total derivative, and for s = 6 there are 
three admissible non-derivative fields. It follows that P2 always satisfy (|58(l . 
and the coefficients in (|59]l can be adjusted to ensure that P4 and Pq satisfy 
(|58(l as well. The role of the boundary condition (J2J) in this counting is in 
reducing the number of independent boundary fields: the boundary condi- 
tions determine the number of redundant boundary fields. As usual, this 
kind of argumentation fails for higher spins, since the number of admissible 
fields grows too fast. Nonetheless, in many other models similar "low-spin 
test" was successful in detecting full integrability, with no known (to us) 
exceptions. Therefore we regard the above argument as a strong indication 
of integrability of the spherical brane model. 

In fact, the coefficients in (j59j) for P4 and Pq can be determined explicitly 
from the commutativity condition ()56|) . While Ii automatically commutes 
with any I[P] of the form (|53j). the condition [13,15] =0 turns out to be 
rather rigid. As it turns, it has only one solution (for the coefficients b, c, d 
in (|59[)) with suitable properties 4 , namely 

63 = (AT + 2)/3, 

6 5 = 3(A + 4)/5, (61) 

c 5 = 7(A + 4)/5, 

d 5 = (N + 4) (36 A + 59)/600 . 

At large the integrals I3 and I5 reduce to their quadratic in X terms, in 

4 In fact, there are exactly three solutions. The other two correspond to previously 
known systems of commuting integrals. One represents the trivial free- field case, where 
all integrals are quadratic in X (obviously, normalization assumed in H59H has to 

be changed to accommodate this case); such integrals can be compatible only with the 
free field boundary conditions (i.e. no nonlinear constraints on Xb(t), and possibly a 
term ~ j (frX|(r) added to the action llUll 'l. Another solution corresponds to the "KdV 
series" , in which the currents P s (it) are composite fields built from the energy-momentum 
tensor T(u) = — 9X • 9X, as was described in 
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agreement with easily established fact that in the N — » oo limit the spherical 
brane model becomes a free theory |38j . Moreover, we have checked that 
with this choice the integrals (|6U|) indeed satisfy (|57f) in the first nontrivial 
order of the 1/N expansion. 

It looks likely that ones I3 is fixed through Eas. (|59|) . (|61|) . the higher-spin 
integrals I2/C-1 with k = 3, 4, 5 . . . are determined uniquely by the commu- 
tativity condition [I s , I3] = 0. Although for higher spins the brute-force 
computation of the commutators becomes difficult, it is possible to deter- 
mine this way the coefficients explicitly written in (j6()|k 



k{k - l){2k + N - 2) 

->2fc-l = 



C2fc-1 



2(2fc-l) 

k{k - l)[k - 2){2k + \){2k + N - 2) 



(62) 



6(2fc- 1) 



7.3 Boundary state 



In view of the structure ()55|) of the space TC, it is convenient to think of 
the boundary state in terms of associated boundary operator. Natural iso- 
morphism between Tp and Tp (the right movers are replaced by the left 
movers) makes it possible to to establish one to one correspondence be- 
tween the states in Tp <8> Tp and operators in Tp . Thus the boundary state 
I B) can be re-interpreted as an operator B : Tp — > Tp. This idea was 
extensively used in the study of conformal boundary conditions since the 
original works j^HHHI!- I n this interpretation the equation (|5H)) reduces to 
the statement of commutativity, 

[B,I s ]=0. (63) 

The boundary state therefore can be written as 

\B)= [ y2B n (P)\n,P)®\nTP), (64) 

where | n , P ) are the (orthonormalized) simultaneous eigenvectors of the 
operators I s in the space Tp, and B n (P) are corresponding eigenvalues of 

B, 

1 \ n, P) =B n (P) \ n, P). (65) 

The eigenvalue -Bo(P) corresponding to the Fock vacuum | , P ) in Tp coin- 
cides with the overlap ( | P ). This structure of | B ) emphasizes importance 
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of the problem of simultaneous diagonalization of the integrals of motion I s . 
Similar problem was addressed in [22] for somewhat simpler model of inte- 
grable boundary interaction - the minimal CFT perturbed by the boundary 
field $1,3. In that case associated system of integrals I s was the KdV series 
(see footnote ^SJl. It was observed in [22J that the eigenvalues of I s are 
related to corresponding eigenvalues of B in the following way. When the 
length of the boundary 2irR — > oo, the operator log IB admits asymptotic 
expansion in terms of the local integrals I s 5 . We expect similar relation to 
hold in the case of the spherical brane model, namely 



1 ~ M IR exp [ - j>2 C k (27rE*) 1 - 2k I 2fc _! | 
^ fc=0 ' 



(66) 



where Mm is the boundary operator associated with the infrared fixed point, 
E* is the energy scale (|14|). I2&-1 with k > 1 are the integrals (|59| ) .I|6U |) . 
I_i = § ^ = R, and the coefficients Ck are just numbers, independent of 
the parameters go, R, and P. The expansion is expected to hold in the 
sense of asymptotic R — > 00 series (note that lafc-l ~ by dimensional 

counting). The dependence of the matrix elements (|66|) on R and P comes 
through the integrals I2&— l- For instance, for the vacuum- vacuum matrix 
element 

/ 2fc _i = (0,P|I 2fc _i|0,P) (67) 
we have, from Eas. (|59l62|) 



R V 24 

i? 3 V 12 320 / v ; 

1 ^6 , ^ + ^4 , (iV + 4)(37jV + 78) ^ 2 



i? 5 V 8 4800 

iV(iV + 4)(143iV + 282) 



+- 



483840 



5 This structure appears in integrable boundary flows down to the "basic" conformal 
boundary, the one which admits no primary boundary fields but the identity. In those 
cases the expansion in I s corresponds to expansion of the infrared effective action in terms 
of descendents of the identity. Clearly, the Dirichlet boundary X_b = is of that kind. 
For integrable boundary flows down to "elevated" fixed points, which admit nontrivial 
boundary primaries, the asymptotic expansions of B at large R should contain "dual" 
nonlocal integrals of motion |41| as well. This is what we would expect to have in the case 
of 9 — tt of the circular brane. 
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and for general k 



where h 2 = — ^j-. In the limit when R — > oo with H = kept fixed, these 
eigenvalues become simply l2k— l ~~ * R (2TrH) 2k , and therefore in (|66|) are 
just the coefficients of the power series expansion of the boundary energy 
S(H) in P|). 



£(#)=£* ]^C fc (^) 2 \ (70) 

fc=0 



At the moment we have only a conjecture about exact form of these coeffi- 
cients in the spherical brane model, 

C k = (-lf (2A) 1 -* ^-^ff (2fc - i)*<i-2A) + A-i ) (71) 

with A = jtzt and k > 1. It is possible to check that this expression agrees 
with the leading order of 1/N expansion in the O(N) spherical brane model. 
The conjecture (|71[) concerns with the coefficients with k > 1, but it does 
not directly apply to the coefficient Co. According to (|70|). the latter relates 
to the boundary energy at zero external field H, S(0) = CqE*. Presently, 
we do not have definite idea what exact value of Co is. Simply setting k = 
in (|71|) is not quite acceptable, since doing so results in complex numbers for 
generic N. On the other hand, it is likely that exact Co is closely related to 
the expression (|71|) with k = (its real part = 2ir cot(-7rA)?). One indication 
is that (|71|) with k = develops a pole at = 2, presumably signifying the 
small-instanton divergence specific to this value of N. Similar pole at N = 3 
is present in the bulk vacuum energy of O(N) sigma model |42j . 
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